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f***^ ■ The thermal relic abundance of species critically depends on the assumed underlying cosmo- 

| logical model. In the case of neutralinos, freeze-out takes place long before Big Bang Nucle- 

pT) ■ osynthesis, which provides the strongest constraint on the evolution of the Hubble parameter 

in the Early Universe. We show that non-standard cosmologies, such as models featuring a 
quintessential scalar field or primordial anisotropies, can lead to large enhancements in the 
neutralino relic abundance, up to six orders of magnitudes. Within these scenarios, super- 
symmetric models with large neutralino annihilation cross sections may account for the whole 
. inferred amount of cold dark matter, yielding on the other hand large indirect detection rates. 
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1 Introduction 



Observational cosmology is accumulating evidences for a picture of the Universe dominated by 
^ ■ a dark sector made up of a non-luminous and non-baryonic component {dark matter) and by a 

gravitationally non-clustering form of energy with negative pressure ( dark energy^. While the 
standard cosmological lore, often dubbed Concordance or ACDM Model nicely fits most of the 
observational features^, particle physics and theoretical cosmology face the compelling questions 
related to the nature and the amount of both dark matter and dark energy. 

The minimal supersymmetric extension of the Standard Model provides a suitable candidate 
for cold dark matter, namely the lightest neutralino. The thermal relic abundance of neutralinos 
£l x h 2 depends both on the particle physics setup, and therefore on the supersymmetric (SUSY) 
model, and on the underlying cosmological model, which affects the freeze-out of species in the 
Early Universe. Models where the role of dark energy is played by a dynamical scalar field 
{Quintessence), or scenarios where extra energy density components are present at the time of 
neutralino freeze-out, may yield significant modifications in the relic abundance of neutralinos^. 
This in turns translates in the intriguing possibility that SUSY models with large neutralino 
annihilation cross sections, and hence with copious dark matter detection rates, give rise to the 
required amount of cold dark matter. 
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2 Cosmological Enhancement of Neutralino Relic Abundance 



The thermal relic abundance of neutralinos is determined by a Boltzmann equation, in the form: 

^ = SHn - {a cS v) (n 2 - n e 2 q ) . (1) 
When the effective annihilation rate of the particles falls below the expansion rate of the Universe, 

n eq {a cS v)(T)<H(T) (2) 

the species ceases to be in thermal equilibrium and freezes-out, locking, after a residual annihila- 
tions transient, its comoving number density, and therefore its relic abundance today. While the 
left hand side (l.h.s.) of eq. (J2J) depends on the details of the particle physics model, the r.h.s. is 
fixed by the cosmological model. The Hubble rate depends in its turn, through the Friedmann 
equation (assuming a spatially flat Universe), on the energy density of radiation, matter and of 
any other component X as 

H 2 (T) = -L- [p r (T) + p m (T) + p x (T)] , (3) 

In the presence of a single spatially-homogeneous scalar field (f) with a potential V (4>), the 
equation of motion of the field is dynamically coupled to eq. ©, and 

In the present study, we take as a representative case the exponential potential^ 

F(0) = M^exp(-A0/M P1 ) . (5) 

The model then depends only on two parameters, namely the exponent A and the Hubble 
rate at the end of inflation, Hi, which fixes the initial ratio of the energy densities of the 
quintessential field and of radiation El. The behavior of p<f>{a) is sketched in fig. ^ (a) for a 
few sample cases at Hi 3> 1: after a period of kination, where p^ ~ a -6 , the scalar field 
undergoes a cosmological constant-like behavior, and then possibly tracks the dominant energy 
density component. For a suitable choice of (\,Hi), p^ may be compatible with Big Bang 
Nucleosynthesis (BBN) constraints and provide the correct amount of dark energy at the present 
epochl**^, while giving a sizeable contribution to the energy density budget at neutralino freeze- 
out (Tf. . ~ m x /20 ~ O(10) GeV). 

To quantify the enhancement, it is convenient to introduce the parameter £ x , defined as 
the ratio £ x = {px/ Pr) {Tf° x ) of the extra-components X and radiation energy densities at 
the freeze-out temperature of the neutralino, for definiteness taken to be that in the absence 
of extra-components. Fig. ^ (b) shows the neutralino relic abundance on a parameter space 
slice of the minimal Anomaly Mediated SUSY Breaking (mAMSB) scenario in the (rn x , 
plane. In mAMSB, where the lightest neutralino x is wino-like, the occurrence of a relic density 
enhancement mechanism is clearly mandatory to allow \ to be the main CDM constituent: 
fig. H (b) shows that, in presence of a dynamical quintessential scalar field, Q x h? may well lie 
within the preferred WMAP range^for any m x . It has been shown^that the relic abundance 
enhancement at £ 3> 1 scales as AQ oc \/£ , and that the maximal quintessential enhancement 
compatible with BBN bounds is AQ = (Q^ — f2 no 0)/^ n o </> ~ 10 6 . 

Other cosmological scenarios where H(T) is significantly larger than in the ACDM model 
at neutralino decoupling range from scalar-tensor b theories^Sl to homogeneous but anisotropic 

'in the case of scalar-tensor theories it has been found that re-annihilations of neutralinos may however reduce 
the expected size of the relic abundance enhancemenP 




Figure 1: (a): Sample cases of Cosmologies with a Quintessence field. The solid lines refer to the sum of the 
radiation and matter energy density components, the other curves to the Quintessence contributions at A = 2.5 
and at three values of Hi. (b): Neutralino relic density enhancement in the minimal Anomaly Mediated SUSY 
Breaking Model, at tan/3 = 50, /i > and mo = 1 TeV and TO3/2 ranging from 32 to 192 TeV, as a function of 
the neutralino mass m x and of the parameter in the quintessential scenario with A = 3.5. The yellow shaded 
region has Q x h 2 > 0.13, while in the green-shaded strip 0.09 < Q x fi 2 < 0.13. Blue dotted lines correspond to 
points at fixed Qh 2 while black dashed lines to points at fixed enhancements AQ. 

cosmologies 15 in this latter framework, for the simplest case of Bianchi-I type space-times, 
where the metric reads 

ds 2 = -dt 2 + R\ {t)[dx x ) 2 + R 2 2 {t){dx 2 ) 2 + R 2 3 (t)(dx 3 ) 2 , (6) 

an effective shear energy density, scaling as p s ~ a~ 6 , rapidly falls off, but may well dominate at 
Tf. .. The initial size of the shear component, which can be quantified through the temperature 
T s at which p s = p r , dictates the value of the parameter £, which again faithfully parameterizes 
the thermal relic abundance enhancement (fig. |2] (a)). 

3 Low Relic Density Models 

Indirect detection signals from neutralino dark matter annihilations products are searched for in 
various channels: neutrinos or muons from \X annihilations in the center of the Sun, of the Earth 
or in the Galactic center, gamma rays from the Galactic center or from other high-density dark 
matter regions, and antimatter (positrons, antiprotons, antideuterons) from the Dark Matter 
Halo. All indirect detection channels linearly depend on the effective neutralino annihilation 
cross section, (a e sv) at T = 0, and inversely on the square of the neutralino mass m x . The first 
quantity is often, though not always, related to the inverse of the neutralino relic abundance, 
which is however sensitive to {<j c rv ) at temperatures around neutralino freeze-out. The naive 
expectation from the Q/1 2 ~ 3 • 10 _27 cm 3 s _1 /(cj e ffv) relation, for the WMAP range, is shown 
in fig. |2 (b). The largest deviations (see the black dots in the figure, outside the naive WMAP 
strip) occur in presence of coannihilations, thresholds or of resonances. Noticeably, the size of 
the mentioned deviations may be as large as more than five orders of magnitude! In the figure 
we also indicate the clusters of points corresponding to a given lightest neutralino type (binos, 
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Figure 2: (a): The relic abundance Q x h 2 and the relative enhancement Af2 x for the mAMSB model at tan/3 = 45, 
fi > 0, mo = 1.5 TeV and m 3( / 2 = 100 TeV as a function of the energy density ratio of non-standard over standard 
components at the temperature of neutralino freeze-out (£). The red lines correspond to quintessential models 
with an exponent A = 3.5, while blue crosses represent models with primordial anisotropies of various initial size, 
(b): A scatter plot of supersymmetric models in the plane (m x , {<T e sv)T=o/'rnt). The various symbols correspond 
to different relic abundance ranges in the standard ACDM cosmological model. The blue line indicates the same 
mAMSB slice as in fig. 0(b), while the black line represents the stau coannihilations strip points at tan/3 = 50, 
Ao — and fj, > in mSUGRA, compatible with the WMAP range. The green band corresponds to the naively 



expected (o a gv) from the rule of thumb Q.h ~ 3 ■ 10 



1 /(ccsv), for Qh 2 in the WMAP preferred range. 



higgsinos, winos). Conspicuous indirect detection signals are generically expected from low relic 
density models, as one can read out from fig. [21(b). Models with wino or dominantly higgsino- 
like lightest neutralinos with m x ~ £7(100) GeV fall within this category (mAMSB, SUGRA 
with non-universal gaugino masses). In this respect, not only does the cosmological thermal 
relic abundance enhancement render low relic density models compatible with the CDM content 
of the Universe, but it also implies promising indirect detection signals at future facilities^. 
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